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Abstract 

Some extensions of the Standard Model provide Dark Matter candidate parti- 
cles with sub-GeV mass. These Light Dark Matter particles have been considered 
for example in Warm Dark Matter scenarios (e.g. the keV scale sterile neutrino, 
axino or gravitino). Moreover MeV scale DM candidates have been proposed 
in supersymmetric models and as source of the 511 keV line from the Galactic 
center. In this paper the possibility of direct detection of a Light Dark Mat- 
ter candidate is investigated considering the inelastic scattering processes on the 
electron or on the nucleus targets. Some theoretical arguments are developed and 
related phenomenological aspects are discussed. Allowed volumes and regions for 
the characteristic phenomenological parameters of the considered scenarios are 
derived from the DAMA/Nal annual modulation data. 
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1 Introduction 

Some extensions of the Standard Model provide Dark Matter (DM) candidate particles 
with sub-GeV mass; in the following these candidates will be indicated as Light Dark 
Matter (LDM). 

Light Dark Matter particles have been considered for example in Warm Dark Mat- 
ter scenarios such as e.g. keV-scale sterile neutrino [1], axino or gravitino [2]. In 
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addition, MeV-scale 1 particles (e.g. axino [4], gravitino [5], heavy neutrinos [6, 7], 
moduli fields from string theories [8], Elko fermions [9]) have been proposed as dark 
matter and as source of 511 keV gamma's from the Galactic center, due either to 
DM annihilation [10, 11] or to decay 2 [4, 7] in the bulge. Moreover, also supersym- 
metric models exist where the LSP naturally has a MeV-scale mass and the other 
phenomcnological properties, required to generate the 511 keV gammas in the galactic 
bulge [13]. 

In this paper the direct detection of LDM candidate particles is investigated consid- 
ering the possible inelastic scattering channels either on the electron or on the nucleus 
target. We note, in fact, that - since the kinetic energy for LDM particles in the galac- 
tic halo does not exceed hundreds eV - the elastic scattering of such LDM particles 
both on electrons and on nuclei yields energy releases well below the energy thresholds 
of the detectors used in the field; this prevents the exploitation of the elastic scatter- 
ing as detection approach for these candidates. Thus, the inelastic process is the only 
possible exploitable one for the direct detection of LDM. 

The following process is, therefore, considered for detection (see Fig.l): the LDM 
candidate (hereafter named vh with mass tuh and 4- momentum fc M ) interacts with 
the ordinary matter target, T, with mass tot and 4-momentum p^. The target T 
can be either an atomic nucleus or an atomic electron depending on the nature of 
the vh particle interaction. As result of the interaction a lighter particle is produced 
(hereafter vl with mass rriL < tuh and 4-momentum k' ) and the target recoils with 
an energy Er, which can be detectable by suitable detectors. 




Figure 1: Inelastic scattering process considered for the detection of the Light Dark 
Matter candidate. The 4-momenta are defined in the laboratory frame. The target T 
can be either an atomic nucleus or an atomic electron. 

The lighter particle vl is neutral and it is required that it interacts very weakly 
with ordinary matter; thus, the vl particle escapes the detector. In particular, the 
vl particle can also be another DM halo component (dominant or subdominant with 
respect to the vh one), or it can simply be a Standard Model particle (e.g. vl can be 
identified with an active neutrino in the scenario of the MeV axino of rcf. [4], where 
the diagrams a + e~ — ► v^ T + e~ have been considered). 

1 lt is worth to note that MeV Dark Matter particles can also be considered for a possible solution 
of the missing satellite problem [3]. 

2 We also remind that the possible decay of DM particles into a lighter state, with MeV scale 
splitting, has been considered in refs. [12]. 
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The production mechanism of the vh (and eventually of the v£) in the early Uni- 
verse are beyond the scope of the present paper and can be investigated in the future; 
however, since mn > tul, it is important to require that the lifetime of the possible 
decay of the vh particle is longer than the age of the Universe. In particular, we note 
that the detection process of Fig.l also offers a possible (unavoidable 3 ) decay channel 
of vjj into i>l and TT pair if > + 2mx- 

Finally, we remind that a LDM particle detection has already been discussed in 
rcf. [14], where the possible light bosonic (keV mass axion-like particles) Dark Matter 
candidate has been studied, considering the total conversion of its mass into electro- 
magnetic radiation. In the present case, because of the presence of the particle in 
the final state, the LDM can be either a boson or a fermion. 

In the following, the direct detection of LDM candidate particles is investigated, 
some theoretical arguments are developed and related phenomenological aspects are 
discussed. In particular, the impact of these DM candidates will also be discussed 
in a phenomenological framework on the basis of the 6.3 a C.L. DAMA/Nal model 
independent evidence for particle Dark Matter in the galactic halo [15, 16]. We remind 
that various corollary analyses, considering some of the many possible astrophysical, 
nuclear and particle Physics scenarios, have been analysed by DAMA itself [14, 15, 16, 
17, 18, 19, 20], while several others are also available in literature, such as e.g. refs. 
[21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. Many other scenarios can be considered as well. 
At present the new second generation DAM A/LIBRA set-up is running at the Gran 
Sasso National Laboratory of the I.N.F.N.. 



2 Detectable energy in the inelastic scattering 

In this section the recoil energy, Er, of the target T for the process of Fig.l is evaluated. 

Considering a Lorentz transformation with (3 — (hereafter fc , k' , pa and 

p' are the time components of the respective 4-momenta in the laboratory frame, see 
Fig.l), it is possible to evaluate the energy conservation relation in the center of mass 
(CM) frame: 

\fs = fc ,CM +P0,CM = k' CM +Pq,CM • (!) 

Knowing that |fc'cM| = Ip'cmI; one obtains: 

P0,CM = m2T 2 jf L = \/p'cM+™T ■ (2) 

Defining the Lorentz boost factor: 7 = l/\/l — (3 2 = ka ^° one can write the total 
energy of the target in the laboratory frame after the scattering by means of a Lorentz 
transformation: 

Po = 7 (po,cm + ■ p'cm) (3) 

Assuming the target at rest before the scattering, i.e. p = and p = rriT, and the 
non-relativistic nature of the LDM particle, i.e. k ~ mjj and k ~ tuhVldm, with 



3 All the other decay channels requires the knowledge/assumptions of some other particle couplings. 
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vldm ~ 10 3 c, the energy released to the target T: Er = ko — k' = p' — po (that is, 
the target recoil energy) can be written as: 

Er = Po-niT = (7Po,cm - m r) + lPp C M cosd CM (4) 
= (E r ) + E+ ~ E - cos9cm (5) 

where: i) 9cm is the scattering angle in the CM frame; ii) (Er) is the average energy; 
iii) E_ is the minimum recoil energy; iv) E + is the maximum recoil energy. 

Considering the non-relativistic nature of the LDM, one gets: s ~ (tot + to#) 2 , 
and for the (Er) average energy: 

, , E + + E- to A 

(Er) = • 6 

There A = tur — tur is the mass splitting of the vr-vr particle system and m = 
mH+m L j g ^ ne avera g e mass f the v H -v L particle system. 

The spread of the recoil energy: 



E +- E - „ k 8WT (7) 

(E R ) y to A(m H + to t ) 

can be appreciable only for hundreds MeV LDM interacting on heavy nucleus targets. 
As an example, assuming uir = 1 MeV and A = 0.5 MeV, the spread in energy is few 
10~ 3 both for electrons and for every nuclear targets. Therefore, since the LDM kinetic 
energy is expected to be negligible with respect to A, one has: (Er) 3> (£+ — E_), 
and the recoil energy is approximatively fixed at the (Er) value. 

In order to offer to the reader just a view of the detectability of LDM particles, 
Fig. 2 shows examples of the configurations in the plane A vs to# (shaded areas), cor- 
responding to released energies in Nal(Tl) within the energy interval 1-6 keV electron 
equivalent. The lower area in Fig. 2 refers to LDM interacting on electron targets, 
while the upper area refers to LDM interacting on the Na and I targets (the effect 
of the quenching factor [15] and of the channeling in the Nal lattice [19] have been 
considered). The dashed line ran = A marks the configurations where vl is a massless 
particle (or also a very light particle, such as e.g. an active neutrino or a nearly mass- 
less sterile one or the light axion, etc.). The configurations characterized by A > 2to 6 
(dark area; hereafter m e is the electron mass) are of interest for the positron anni- 
hilation line from the galactic center due to possible (energetically allowed) decay: 
vr — > VLe + e~ . Finally, the configurations on the right of the vertical line at to# = 
511 keV are of interest for the possible annihilation process: vrvr — > e + e~ in the 
galactic center, while the configurations of to# on the right of the solid or dotted 
curves are of interest for the possible annihilation processes: vrDl — » e+e~ (solid 
curve), v^vr — > e+e~ (solid curve), and vl&l — * e + e~ (dotted curve) in the galactic 
center. 
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Figure 2: Configurations in the plane A vs ran (shaded areas), corresponding to 
released energies in Nal(Tl) within the energy interval 1-6 keV electron equivalent. 
The upper area is due to LDM interaction on Na and I nuclear targets, while the lower 
area to LDM interaction on electron target. The dashed line (ran — A) marks the case 
where vl is a massless particle. The configurations characterized by A > 2m e (dark 
area) are of interest for the positron annihilation line from the galactic center through 
the decay: vh — > v^e^er . The thresholds of the possible annihilation processes: 
vrvh — > e + e~ (solid vertical line at m# = 511 keV); vh&l — * e + e~ (solid curve); 
vlVh —* e + e~ (solid curve); vl&l — * e + e~ (dotted curve) are shown. 

3 Interaction rate of LDM particle 

The interaction rate of the LDM particle with the target T and for the studied process 
can be written as: 



where: i) p VH — £po, with po local halo density and £ < 1 fractional amount of vh 
density in the halo; ii) T]t is the target number density in the detector; iii) is 
the differential cross section of the considered LDM particle inelastic scattering with 
the target T; iv) f(v) is the velocity (v) distribution of the vjj particles in the Earth 
(laboratory) frame. 

Since the sub-GeV LDM wavelength (A = § > 10 3 fm) is much larger than the 
nucleus size, the targets can be considered as point-like and the form factors of the 
targets can be approximated by unity. In such a reasonable hypothesis and assuming, 




(8) 
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for simplicity, the isotropy of the differential cross section , one gets: 



dE R — E_ 



-Q{E+ - E R )e(E R - EJ) , (9) 



where: i) <jt is the target cross section; ii) the Heaviside functions, 9, define the 
domain of the differential cross section. We note that E + and E_ depend on the 
target mass and slightly on the LDM velocity, v. 

Generally, a? can be a function of v, depending on the peculiarity of the particle 
interaction. In the following, we adopt a widely considered approximation for the 
non-relativistic case [11, 13, 31] (the subscript T drops): 

<jv ~ a + bv 2 . (10) 

There a and b are constant and depend on the peculiarity of the particle interaction 
with the target T; they will be considered as free parameters in the description of the 
process. In the following, in order to deal the direct detection process with more usual 
parameters, the cross sections <Jq = and a^ n = bv will be used as free parameters; 
they are related to the a and b parameters rescaled with the Dark Matter local velocity, 
v a , [32, 33]. 

Thus, eq. (8), becomes: 

dRr Pu H f, T 2 , t 2^ &{E+-E R )@{E R -E_) 



dE R vomH 



f t T 2 . T 2\ _ E R )Q(E R - E-) 3 

J \ a o v o + a m v ) E+-E- f( v ) dv - ( n ) 



Since: J °° e ( g + Er)&{e r e _ the total rate Rtf* can be written as: 

Rt ot = r dRr dER = Vt J^_ {a r vl + a r {v2)) (12) 
Jo dtj R v rriH 

Similarly as the light bosonic case discussed in ref. [14], the a^(v 2 } term provides 
an annual modulation of the expected counting rate for LDM interactions. In fact, 
the velocity of the LDM particle in the galactic frame can be defined as: v g = v + v®. 
Thus, one obtains for non rotating halo: (v 2 ) = (v 2 ) + v^. Since the Earth velocity in 
the galactic frame (v®) is given by the sum of the Sun velocity (v Q ) and of the Earth's 
orbital time-dependent velocity around the Sun (vsE(t)), neglecting the v 2 SE term one 
gets: 

(i; 2 ) ~ (v 2 ) +v 2 Q + 2v Q ■ v SE (t) ~ (v 2 ) + v% + v G v SE cos{uj(t - to)) ■ (13) 

In the last equation the angle of the terrestrial orbit with respect to the galactic plane 
has been considered to be ~ 60°. Moreover, we consider u> — 2ir/T with T = 1 
year, and the phase to (which corresponds to ~ 2nd June, that is, when the Earth's 
speed in the galactic frame is at the maximum). The Sun velocity can be written as 
| if© | ~ Vo + 12 km/s, where i>o is the local velocity, whose value is in the range 170-270 
km/s [32, 33]. The Earth's orbital velocity is v$e — 30 km/s. Finally, (v 2 ) depends 



on the halo model and on the vo value (just for reference, in the particular simplified 



case of isothermal halo model: (v 2 ) = ffp)- 



4 Different assumptions would produce quite similar phcnomcnologies, since — as mentioned 
(En) 2> — E-), and the recoil energy is approximatively fixed at the (En) value. 
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In conclusion, the expected signal is given by the sum of two contributions: one in- 
dependent on the time and the other one dependent on the time through cos(uj(t — to)), 
with relative amplitude depending on the adopted scenario. 



3.1 Interaction with atomic electrons 

For atomic electron targets, the energy interval (E_,E + ) is always very narrow 
^fci < /em %j . Therefore, using the approximation e ( g +~-^«K-X g H--E-) ^ 
S(E R — (En)), the interaction rate with atomic electrons can be simply written as: 

§L „ Ve J^_ ( a e v 2 + a e m{v 2 }) ^ {Er)) (u) 

where r\ e is the electron number density in the target detector and ct§ and are the 
average (over all the atomic electrons) cross section parameters for electrons in the 
target material. 

After the interaction the final state can have - beyond the vl particle - either a 
prompt electron and an ionized atom or an excited atom plus possible X-rays/ Auger 
electrons. Therefore, the process produces X-rays and electrons of relatively low en- 
ergy, which are mostly contained with efficiency ~ 1 in a detector of a suitable size. 

Therefore, the differential counting rate can simply be written accounting for the 
detector energy resolution, by means of the G(E,E R ) kernel, which generally has a 
Gaussian behaviour: 



'!^ = Jg(e,e, — 7 = /?,"" x — , — ^— . do) 



dR e ,th root 1 ( E - <E "» 2 



There S is the energy resolution of the detector and is the area of the Gaussian 
peak centered at the (E R ) value. 

Finally, the expected differential rate is given by the sum of two contributions: 
^jgf = So + S m -cosuj(t — to), where So and S m are the unmodulated and the modulated 
part of the expected differential counting rate, respectively. 



3.2 Interaction with nuclei 

As regards the interaction with target nuclei, when considering the cases of heavy 
target and for relatively "high" LDM mass the energy interval (E^,E + ) cannot be 
neglected. Therefore, the interaction rate with the nucleus T can be approximated as: 

dRr Pu H / t 2 I t, 2 6(El-E R )e(E R -ET) 

~dE R * ^ V ° + ^ } ) Ej^W ' (16) 

where t)t is the nucleus number density in the material and E± are the E± values for 
the nucleus T calculated at v = (v). 

Considering the case of a crystal detector target (as e.g. the Nal(Tl)), the possible 
channeling effect [19] has also to be taken into account. The detection of a nuclear recoil 
of kinetic energy, E R , is, in particular, related to the response function (Edet\E R ) 
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[19], where Edet is the released energy in keV electron equivalent. Thus, the differential 
Edet distribution for the nucleus T when accounting for the channeling effect can be 
written as: 




(E de t) = J ^-(E det \E R )^(E R )dE 



(17) 



dEdet 



R ■ 



Hence, the expected differential counting rate for LDM interaction with nucleus T 
can easily be derived accounting for the detector energy resolution kernel G(E, Edet)' 



Finally as for the electron interacting LDM, the expected differential rate is given 
by the sum of two contributions: = So + S m ■ cosu(t — to), where So and S m are 
the unmodulated and the modulated part of the expected differential counting rate, 
respectively. 

3.3 Some examples on interaction rate in Nal(Tl) 

In this section as template purpose, some examples on interaction rate of LDM particles 
in Nal(Tl) detectors are given. In particular, since the expected energy distribution 
for electron 5 interacting LDM has a simple Gaussian behaviour (sec cq. (15)), only 
some examples for the nucleus interacting LDM will be shown in the following. 

In Nal(Tl) detectors, both Na and I nuclei {rjN a — r\i ~ 4 x 10 24 kg -1 ) can 
contribute to the detection of LDM. Therefore, the total expected counting rate is given 
by the sum of two contributions similar as eq. (18), arising from the two nuclei species 
(T = Na, I). The cross section parameters Oq and tr^ for Na and I are related with 
some model dependent scaling laws; in particular, in the following, the two illustrative 
cases of coherent (a coh ~ 4^— ~ 4^) and incoherent (p mc ~ a Na ~ a 1 ) nuclear scaling 
laws are investigated. Obviously, other different nuclear scaling laws are in principle 



As template purpose, in Fig. 3, the expected energy distributions of the unmod- 
ulated/modulated part of the differential rate in Nal(Tl) are shown for interactions 
of LDM with mjj = 100 MeV on the target nuclei. In particular, there are shown 
the cases of coherent cross section scaling laws (Fig. 3 -left) with A = 4.8 MeV and 
£ct§ o/i < £o% h ~ 2 x 10~ 6 pb and the case of incoherent cross section scaling laws 
(Fig. 3 -right) with A = 0.95 MeV and £<r£ nc < £<r™ c ~ 20 x 10~ 3 pb. The A5 halo 
model (a NFW halo model with local velocity equal to 220 km/s and density equal 
to the maximum value, see rcf. [15, 16]) has been considered. The quenching factors 
have been assumed as the case A of ref. [15, 16] and the channeling effect has been 
included, as mentioned. 

It is worthwhile to note that similar behaviors can also be obtained by using other 
choices of the halo model, quenching factors and values of the masses of the involved 
particles in the interaction. 

5 Here, we only remind that r/ e ~ 2.6 X 10 26 kg^ 1 in Nal(Tl). 




(18) 



possible. 
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Figure 3: Examples of the energy distributions of the unmodulated (solid) and of the 
modulated (dotted) parts of the expected differential rate in Nal(Tl) for interactions 
of LDM with ran = 100 MeV on the target nuclei. Left: case of coherent cross section 
scaling laws; here A = 4.8 MeV and £ct§ o/i < £a% h ~ 2 x 10~ 6 pb. Right: case 
of incoherent cross section scaling laws; here A = 0.95 MeV and £a™ c <C C' 7 ™ — 
20 x 10~ 3 pb. The A5 halo model (a NFW halo model with local velocity equal to 220 
km/s and density equal to the maximum value, see ref. [15, 16]) has been considered. 
The quenching factors have been assumed as the case A of rcf. [15, 16]. The channeling 
effect has been included; see text. The vertical dotted lines correspond to the energy 
threshold of the Nal(Tl) detectors used in DAMA/Nal set-up. 

4 Data analysis and results for LDM candidates in 
DAMA/Nal 

The 6.3 a C.L. model independent evidence for Dark Matter particles in the galactic 
halo achieved over seven annual cycles by DAMA/Nal [15, 16] (total exposure ~ 1.1 x 
10 5 kg x days) can also be investigated for the case of a LDM candidate (in addition 
to the other corollary quests analyzed by DAMA itself [14, 15, 16, 17, 18, 19, 20] and 
available in literature, e.g. refs. [21, 22, 23, 24, 25, 26, 27, 28, 29, 30]). 

In the following, the same dark halo models and related parameters given in table 
VI of ref. [15] are considered. The related DM density is given in table VII of the 
same reference. Note that, although a large number of self-consistent galactic halo 
models have been considered, still many other possibilities exist. As regards the case 
of LDM interacting on nuclei, also the Na and I quenching factors [15] and the Na and 
I channeling response functions given in [19] have been considered. The uncertainties 
on the Na and I quenching factors have been taken into account as done in ref. [15]. 
In addition, the presence of the existing Migdal effect and of the possible SagDEG 
contributions - we discussed in [18, 17], respectively - will be not included here for 
simplicity. 

As regards the cross section parameters, the a terms do not contribute to the 
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annual modulation of the signal. In particular, the measured modulation amplitude 
easily allows to derive <7q < <J m - In the following, for simplicity, the contributions of 
the a m terms are assumed to be dominant with respect to the contributions of the ao 
ones. 

The results are calculated by taking into account the time and energy behaviours of 
the single-hifi experimental data through the standard maximum likelihood method 7 . 
In particular, they are presented in terms of slices of the three-dimensional allowed 
volume (m H , A, £,<r m ) - where a m is alternatively a^, a^ 1 and a™ c - obtained as 
superposition of the configurations corresponding to likelihood function values distant 
more than 4cr from the null hypothesis (absence of modulation) in each one of the 
several (but still a very limited number) of the considered model frameworks. In 
this way one accounts for at least some of the existing theoretical and experimental 
uncertainties. It is worth to note that the inclusion of other existing uncertainties 
would further extend the allowed volumes/regions and increase the sets of obtained 
best fit values. 

The projection of the whole 4<r allowed volume on the plane (m^r, A) gives typ- 
ical patterns similar as those reported in Fig. 2. Moreover, since the two regions are 
disconnected, the LDM detection is always dominated by only one of the different 
target contributions. Therefore, in the following, slices of the allowed volume either 
for electron interacting LDM or nucleus interacting LDM are presented separately. 

4.1 Case of electron interacting LDM 

In case of electron interacting LDM, the projection of the 4er allowed volume on the 
plane (rriH, A) for the same dark halo models and parameters described in ref. [15] 
is reported in Fig. 4. The allowed vhh values and the splitting A are in the inter- 
vals 35 keV < m H < O(GcV) 8 and 1 kcV < A < 80 kcV, respectively. It is worth to 
note that in such a case the decay through the detection channel: uh ~^ vi,e + e~ , 
is energetically forbidden for the given A range. The configurations with mu > 511 
keV (dark area in Fig. 4) are instead of interest for the possible annihilation processes: 
v rVh — * e + e~, vh^l — * e + e~, vl^h — > e + e _ , and vl^l — > e + e~, in the galactic 
center. 

As examples, some slices of the 3-dimcnsional allowed volume for various ttih values 
in the (£o^ vs A) plane are depicted in Fig.5-ie/t in the scenario given above. 

6 The single-hit events are those events where only one detector of many actually fires in a multi- 
detectors set-up. In particular, the vn particle is no more present after the inelastic interaction and, 
therefore, can only provide single-hit events in a multi-dctcctors set-up. 

7 Shortly, the likelihood function is: L = Iljjj.e Mi J fc ' jfc , , where Nijk is the number of events 

ijk • 

collected in the i-th time interval, by the j-th detector and in the fc-th energy bin. follows a 

Poissonian distribution with expectation value fiijk = [bj k + Sq ^ + S m ^ ■ cosw(ti — to)]Mj AtjAEe^j.. 
The unmodulated and modulated parts of the signal, So k and S m ucosu)(ti — to), respectively, are 
here functions of the LDM mass van, the splitting A and the cross section <r m . The bj^ is the 
background contribution; Atj is the detector running time during the i-th time interval; AE is the 
energy bin; ejk is the overall efficiency and Mj is the detector mass. 

8 For values of ran greater than O(GeV), the definition of Light Dark Matter is no more appro- 
priate. Moreover, the kinetic energy of the particle would be enough for the detection also through 
the elastic scattering process, as demonstrated in ref. [20]. 
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Figure 4: Case of electron interacting LDM. Projection of the 4tr allowed 3-dimensional 
volume on the plane (run , A) for the same dark halo models and parameters described 
in ref. [15];see text. The dashed line (mjj = A) marks the case where vl is a massless 
particle. The decay through the detection channel, vjj — * ^Le + e~, is energetically not 
allowed for the selected configurations. The configurations with mu > m e (dark area) 
are interesting for the possible annihilation processes: vh^h — * e + e _ , vjjvl — ► e + e~, 
vl&h — > e + e~, and Vl^l - > e + e~ in the galactic center. The three nearly vertical 
curves are the thresholds of these latter processes as mentioned in the caption of Fig. 2 
and in the text. 




The slice of the 4<r allowed 3-dimensional volume for run = A is shown in Fig. 5- 
right. This slice has been taken along the dotted line of Fig. 4, restricting mi ~ 0, that 
is for a massless or a very light vl particle, such as e.g. either an active neutrino or a 
nearly massless sterile one or the light axion, etc. 



Finally, it is worth to summarize that electron interacting LDM candidates in the 
few-tens-keV/sub-MeV range are allowed (see Figs. 4 and 5). This can be of interest 
for example in the models of Warm Dark Matter particles - such as e.g. Weakly 
Sterile Neutrino [1], Axino or Gravitino [2] - or in the models where the DM is made 
of Moduli fields [8]. Moreover, also configurations with mjj in the MeV/sub-GeV 
range are allowed; similar LDM candidates (such as e.g. axino [4], sterile neutrino 
[7] , moduli fields from string theories [8] , and even MeV-scale LSP in supersymmetric 
theories [13]) can also be of interest for the production mechanism of the 511 keV 
gammas from the galactic bulge [10, 11]. 
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Figure 5: Case of electron interacting LDM. Left: examples of some slices of the 4<r 
allowed 3-dimensional volume for various mu depicted in the (£(7^ vs A) plane. Right: 
slice of the 4tr allowed 3-dimensional volume for mu = A, that is for a massless or a 
very light vj, particle, as e.g. either an active neutrino or a nearly massless sterile one 
or the light axion, etc. The same dark halo models and parameters described in ref. 
[15] have been used; see text. 



4.2 Case of nucleus interacting LDM 

In case of nucleus interacting LDM, the projections of the 4er allowed volumes on the 
plane (mjj, A) are reported in Fig. 6 for the two above-mentioned illustrative cases 
of coherent (left) and incoherent (right) nuclear scaling laws. They are evaluated for 
the same dark halo models and parameters described in ref. [15, 19]. The allowed 
ran values and the splitting A are in the intervals 8 MeV < run < O(GeV) 9 and 
25 keV < A < 120 MeV, respectively (see Fig. 6). It is worth to note that in such a 
case the decays through the diagram involved in the detection channel (e.g. in nucleon 
anti-nucleon pairs or in meson(s), as vh — > vl^°) are obviously energetically forbidden. 
Moreover, there are allowed configurations that could contribute - in principle, if 
suitable couplings exist - to the positron generation in the galactic center; in fact, the 
decay vu —* v^e^e^ is energetically allowed for A > 2m e (dark area in Fig. 6), while 
the annihilation processes into e + e~ pairs are energetically allowed for almost all the 

9 We remind that for mu values greater than O(GeV) the detection would also be possible through 
the elastic scattering process [15, 16, 17, 18, 19]. 
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Figure 6: Case of nucleus interacting LDM. Projections of the 4tr allowed 3-dimensional 
volumes on the plane (run, A) for the two above-mentioned illustrative cases of co- 
herent (left) and incoherent (right) nuclear scaling laws. They are evaluated for the 
same dark halo models and parameters described in ref. [15, 19]. The regions enclose 
configurations corresponding to likelihood function values distant more than 4cr from 
the null hypothesis (absence of modulation). The dashed lines (mjj = A) mark the 
case where i>l is a massless particle. The decays through the diagram involved in the 
detection channel are energetically forbidden. See text. 



allowed configurations. 

It is worth to note that for nuclear interacting LDM the 3-dimensional allowed con- 
figurations are contained in two disconnected volumes, as seen e.g. in their projections 
in Fig. 6. The one at larger A at mu fixed is mostly due to interaction on Iodine 
target, while the other one is mostly due to interaction on Sodium target. 

As examples, some slices of the 3-dimensional allowed volumes for various mn 
values in the (£c m vs A) plane are depicted in Fig. 7 for the two above-mentioned 
illustrative cases of coherent (left panel) and incoherent (right panel) nuclear scaling 
laws. 

The slices of the 4cr allowed 3-dimensional volumes for m# = A are shown in 
Fig. 8 for the two illustrative cases of coherent (left panel) and incoherent (right panel) 
nuclear scaling laws. These slices have been taken along the dotted lines of Fig. 6, 
restricting tol — 0, that is for a massless or a very light vl particle, as e.g. either an 
active neutrino or a nearly massless sterile one or the light axion, etc. 

Finally, it is worth to summarize that LDM candidates in the MeV/sub-GeV range 
are allowed (see Figs. 6, 7 and 8). Also these candidates, such as e.g. axino [4], 
sterile neutrino [7], moduli fields from string theories [8], and even MeV-scale LSP 
in supersymmetric theories [13], can be of interest for the positron production in the 
galactic center [10, 11]. 
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Figure 7: Case of nucleus interacting LDM. Examples of some slices of the 3- 
dimensional allowed volumes for various m# values in the (£<r m vs A) plane for the 
two above-mentioned illustrative cases of coherent (left) and incoherent (right) nu- 
clear scaling laws. The 3-dimensional volumes enclose configurations corresponding to 
likelihood function values distant more than 4<r from the null hypothesis (absence of 
modulation). The same dark halo models and parameters described in ref. [15, 19] 
have been used. 
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Figure 8: Case of nucleus interacting LDM. Slices of the 4cr allowed 3-dimensional 
volumes for mn = A, that is for a massless or a very light vl particle, as e.g. either 
an active neutrino or a nearly massless sterile one or the light axion, etc. They are 
evaluated for the two illustrative cases of coherent (left panel) and incoherent (right 
panel) nuclear scaling laws, using the same dark halo models and parameters described 
in ref. [15, 19]. 
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5 Conclusions 



In this paper, the possibility of direct detection of a Light Dark Matter candidate 
particle has been investigated. The inelastic scattering processes on the electron or on 
the nucleus targets have been considered; these are the only possible processes useful 
for the direct detection of such LDM candidates. 

Some theoretical arguments have been developed and related phenomenological 
aspects have been discussed. In particular, the impact of the LDM candidate has also 
been analyzed in a phenomenological framework on the basis of DAMA/Nal annual 
modulation data. Allowed volumes and regions for the characteristic phenomenological 
parameters have been derived in the considered model. The allowed phenomenological 
parameters (masses and cross sections) of LDM can be of interest for various LDM 
candidates proposed in theories beyond the Standard Model and for the production 
mechanism of the 511 keV gammas from the galactic bulge. 

In conclusion this paper has shown that - in addition to other candidates, as 
WIMP/WIMP-like particles and axion-like bosons, already discussed by DAMA col- 
laboration elsewhere [14, 15, 16, 17, 18, 19, 20] - there is also possibility for a LDM 
candidate interacting either with the electrons or with the nuclei to account for the 
6.3(7 C.L. model independent evidence for the presence of a particle DM component 
in the galactic halo. 
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